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ABSTRACT 

We have carried out a survey for 'giant pulses' in 6 young, Vela-like pulsars. In no 
cases did we find single pulses with flux densities more than 10 times the mean flux 
density. However, in PSR B1706-44 we have detected giant micro-pulses very similar 
to those seen in the Vela pulsar. In PSR B1706-44 these giant micro-pulses appear on 
the trailing edge of the profile and have an intrinsic width of ~1 ms. The cumulative 
probability distribution of their intensities is best described by a power-law. If the 
power-law continues to higher intensities, then 3.7 x 10® rotations are required to 
obtain a pulse with 20 x the mean pulse flux. This number is similar to the giant pulse 
rate in PSR B1937+21 and PSR B1821-24 but significantly higher than that for the 
Crab. 
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1 INTRODUCTION 

Despite three decades of intensive study, the mechanism pro- 
ducing pulsar radio emission is poorly understood. Fluctu- 
ations in the intensity of the radio radiation provide im- 
portant constraints on plausible mechanisms. Single-pulse 
studies of bright pulsars detect a variety of patterns in the 
intrinsic intensity fluctuations, including nulling and drift- 
ing phenomena. The distribution of integrated pulse ener- 
gies, however, has only a modest dispersion. Johnston et al. 
(2001) showed that in the Vela pulsar, 99.5% of all pulses are 
within a factor of 3 of the mean flux density, {S) and that 
the histogram of pulse energies is a Gaussian when plotted 
i n the log. This distribution seems typical for most pulsars 
(|Hesse fc Wielebinski 1974]; [Ritchings 197^). 



In contrast, the Crab pulsar emits pulses wit h flux den- 
sities > 20 X (S), extending up to > 2 x 10^(5') (Lundgren 



et al. 1995) which were instrumental in the o riginal detec- 
tion of the Crab ( Staelin fc Reifenstein 1968 ) . These giant 
pulses are typically broadband ( |MofictrT997 ; Ballmen et al 



199£ ) and of short duration, with widths of order a few /is 



and structure down to 10 ns (Hankins 1996). They are local 



ized to the main and interpulse phase windows and follow 
an intensity distribution best characterized as a power law 
with index ~ 3 — 3.5. 

The discovery of si milar pulses from the m i llisecond 
pulsar PSR B1937+21 (pallmen fc Backer 1995|; pognard 



et al. 1996) was surprising. The pulses are extremely short 
(t < 0.3/is) events conflned to small phase windows trail- 
ing the main pulse and interpulse, again with a n approxl 



known radio pulsar with an estimated magnetic field at the 
light cylinder larger than that of the Crab, it has been sug- 
gested that this is a key parameter co ntrolling the giant 
pulse phenomenon ( Cognard et al. 1996| ). 



mately power-law distr ibution of pulse energies (Kinkhab 
wala fc Thorsett 2000|). Since PSR B1937+21 is the only 



Johnston et al. (2001) have recently found that a small 
subset of pulse phases from the Vela pulsar have a very wide 
distribution of peak fluxes. The pulses are localized to a 
phase window ~1 ms prior to the bulk of the integrated 
pulse emission, are of short duration and are highly polar- 
ized. Johnston et al. (2001) called these giant micro-pulses. 
It is not clear if these events are related to true giant pulses; 
the largest Vela pulses observed to date have S < 10{S), but 
these narrow pulses have peak fiuxes exceeding 40 x the in- 
tegrated peak intensity. Kramer, Johnston & Van Straten 
(2001) have shown that these giant micro-pulses have a 
power-law distribution and the extended tail of the distribu- 
tion may continue into the true giant pulse regime. Cairns, 
Johnston & Das (2001), in contrast, recently showed that a 
log-normal distribution provided an excellent fit to the fiux 
densities in individual phase bins across the main peak of the 
Vela profile. It seems likely that the same distribution is also 
applicable to other pulsars. Therefore a potential discrimi- 
nator of giant pulse activity is the change from a log-normal 
distribution to a power-law one. 

To explore the connection between the giant pulses in 
the Crab and PSR B1937+21 and the large individual pulses 
in the Vela pulsar, we obtained fast time-sampled data for 
s everal young and millisec ond pulsars. In an earlier paper 
(Romani & Johnston 2001) we reported the detection of gi- 
ant pulses from the millisecond pulsar with the next highest 
known light cylinder field, PSR B1821-24. In this paper 
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Table 1. Pulsars observed in the survey 



Name 




P 


DM 


Age 


Tobs 


Npuls 






(ms) 


(cm~'^pc) 


(kyr) 


(hr) 




B1046- 


-58 


123.7 


129.1 


20.3 


3 


81785 


J1105- 


6107 


63.2 


270.6 


63.4 


3 


168099 


J1420- 


6048 


68.2 


359.8 


13.0 


4 


205222 


B1509- 


-58 


150.7 


253.2 


1.6 


3 


61609 


J1617- 


5055 


69.4 


467.0 


8.0 


3 


153536 


B1706- 


-44 


102.4 


75.7 


17.5 


3 


95536 



we report on our survey of young pulsars focussing on the 
results obtained for PSRs B1046-58 and B1706-44. 
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Figure 1. Two large single pulse from PSR B1046— 58 shown 
together with the integrated profile magnified by a factor of 10. 



2 OBSERVATIONS AND DATA ANALYSIS 

Table 1 lists the relevant parameters for the pulsars ob- 
served in this survey. The pulsars were chosen on the cri- 
teria of small spin periods and/or small age coupled with 
a reasonably high flux density at 1500 MHz. The link be- 
twee n high energy emission an d giant pulses seems a strong 
one (Romani & Johnston 2001), hence we also chose pulsars 



which have been detected at high energies. PSRs B1046-58, 
B1706-44 and J1105-6107 are not known to pulse in the 



rounding nebula ( Pivovaroff , Kaspi & Gotthelf 2000 ; 


Becker, 


Brazier & Triimpcr 1995t Gotthelf & Kaspi 199^). The first 


two are are pulsed in 7-rays (|Kaspi et al. 200C ; 


Thompson 



et al. 1992), the latter is a 7-ray de tection but a sear ch for 



pulses has as yet proved unsuccesful ( Kaspi et al. 2000 ) . PSR 



B 1509-58 is pulsed at both X-ray and 7-ray wavelengths 



(Ulmer et al. 1993 



Rots et al. 



199^ ). 



PSRs J1420-6048 and 



J1617-5055 are pulsed ii i the X-ray and may also be associ 
ated with 7-ray sources ([Roberts, Romani fc Johnston 2001 



Gotthelf, Petre fc Hwang 1997D 

Observations were made on 20-22 May 2001, with the 
Parkes 64-m radio telescope. We used the center beam of 
the 21-cm multi-beam system at an observing frequency of 
1517.5 MHz. The receiver has a system equivalent flux den- 
sity of 30 Jy on cold sky. The back-end consisted of a filter- 
bank system containing 512 channels per polarization each 
of width 0.5 MHz for a total bandwidth of 256 MHz. The 
polarization pairs are summed, each output is then sampled 
at 80 /xs, one-bit digitised, and written to DLT for off-line 
analysis. 

The data were then de-dispersed at the pulsar's nomi- 
nal dispersion measure (DM). The mean and rms of groups 
of 8192 samples (0.65 s) were examined and those which 
showed obvious signs of interference were discarded. The 
data could then be folded synchronously with the pulsar's 
topocentric period to produce a pulse profile. Our nominal 5- 
a sensitivity in 80 /xs is 1.3 Jy. In practice, even after removal 
of high sigma points (clipping), we experienced substantially 
larger background fluctations and the rms exceeded our ex- 
pected rms by a factor of 1.5-3. 

To confirm that we could detect conventional giant 
pulses, we made short observations of the Crab pulsar, dur- 
ing which giant pulses were detected with high significance. 
Short integrations with the Vela pulsar off-axis, producing 



an effective contirmum fiux of only 0.6 mjy, confirmed that 
we could detect the largest individual Vela-type giant micro- 
pulses at faint continuum flux levels. Finally, we observed 
PSR B1937-I-21, obtaining 1784 s of integration (1.15 x 10*' 
pulses). Our sampling provides only 19 bins across the pulse 
profile, but we clearly detect the giant pulse distributions in 
both the main and interpulse components. In both compo- 
nents the giant pulse peaks occur ~ 1 sample after the cor- 
responding peak of the integrated pulse profile. The largest 
main component giant pulse obtained had an energy 445 
Jy/xs, the third largest had 230 Jy/iS. The phasing and in- 
tensity distribution are well matched to the 1.4 GHz results 
reported by Kinkhabwala & Thorsett (2000). 



3 RESULTS 

3.1 PSR B1046-58 

PSR B1046-58 has a mean flux density of 5.9 mJy. It has a 
narrow pulse profile with a peak fiux of 150 mJy, the highest 
of any of the pulsars in our sample. We re-sampled the 80 ^s 
data to produce 256 phase bins across the pulsar period and 
examined the ~82000 single pulses. The rms per phase bin is 
85 mJy; in ~82000 pulses we expect the largest noise sample 
to be 4.3(T. 

On both the leading edge and trailing edges of the pulse 
we see a number of large amplitude pulses, where the peak 
flux exceeds 20 x the integrated flux in these phase bins. Ex- 
amples of these pulses are shown in Figure 1. This behaviour 
is reminiscent of the Vela pulsar, where Krishnamohan & 
Downs (1983) showed that the strong pulses were confined 
to the rising edge of the profile. Kramer et al. (2001) showed 
that the distribution of fiuxes on the rising edge of the Vela 
pulsar has a log-normal distribution with a relatively large 
sigma of 0.6 in the log, in contrast to the giant micro-pulses 
which occur even earlier in phase and which have a power- 
law distribution of fluxes. In PSR B1046-58 the situation 
is less extreme. Figure 2 shows the distribution of fluxes on 
the rising edge of the pulse and their cumulative probability 
distribution. These are consistent with a log-normal fit with 
a sigma of 0.33 in the log. The trailing edge of the pulse 
shows rather similar behaviour. The best fit to these data 
is also log-normal with slightly larger sigma of 0.45 in the 
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Figure 2. Top panel: Histogram of flux densities on tlie leading 
edge of the pulse of PSR B1046-58. The solid line denotes the 
best fit to the data after convolving a gaussian noise distribution 
of width 32.0 mjy with a log-normal pulsar distribution with 
mean of 36.5 mJy and sigma (in the log) of 0.33. Bottom panel: 
Cumulative probability distribution of the data (binned) and best 
fit (solid line). 



log. In contrast, the centre of the pulse is best fit with a 
log-normal distribution with sigma of 0.25. 

The key issue here is that log-normal distribution pro- 
vides an excellent fit to the data at all pulse phases (Cairns 
et al. 2001). There is no evidence of any power-law distribu- 
tion at any phase. The variations seen in sigma across pulse 
phase is consistent with the statement that the modulation 
index is large towards the outside of the pulse and small 
in the centre, a trend which has been known since the mid- 



1970s (Taylor, Manchester & Huguenin 1975; Krishnamohan 



Dow IS 198;) 



3.2 PSR J1105-6107 

The integrated fiux density of this pulsar is 1.49 mJy. The 
pulse profile is double-peaked with a peak fiux density of 60 
mJy. We subdivided the pulse period into 256 phase bins 
for a time resolution of 250 fis and examined the ~ 168000 
pulses. The effective rms was 100 mJy per phase bin. Only 
phase bins towards the centre of the pulse had a flux density 
in excess of that expected from the noise fluctuations. The 
distribution of intensities in the peak flux bin was consistent 
with a log-normal distribution with sigma of 0.3 in the log 
convolved with gaussian statistics for the noise. There was 
no evidence of a power-law tail to the distribution in any 
phase bin. 




Pulse Phase 



Figure 3. A large single pulse from PSR B1706-44 shown to- 
gether with the integrated profile magnified by a factor of 10. 
Virtually all the large pulses are located at this pulse phase. 



3.3 PSR J1420-6048 

The pulsar was observed for a total of 4 hr, split into two 
separate observations. The integrated fiux density for both 
observations is only 0.95 mJy. The pulsar has large duty 
cycle; as a result its peak flux is only 12 mJy. No large sam- 
ples were detected at any pulse phase. Subdividing the pulse 
into 256 phase bins and examining the single pulses yielded 
statistics which were consistent with the noise level of 90 
mJy. One pulse out of ~210000 had, however, a peak flux of 
800 mJy, this is the only sample in the entire data set with 
a flux greater than 5 times the noise rms (450 mJy). This 
high point is located on the leading edge of the second peak, 
where the integrated flux density is only 2.1 mJy. However, 
re-examination of individual 80 /is samples, revealed that 
this high flux is entirely due to one sample. As the DM of 
the pulsar ensures that any real signal should be smeared by 
~7 samples, we therefore believe this does not originate from 
the pulsar but is likely to be interference or an instrumental 
glitch. 



3.4 PSR B1509-58 

The integrated flux density during our 3 hr observation was 
only 0.96 mJy. This pulsar has a broad single component 
with a peak fiux density of only 8 mJy. The DM is large, 
effectively smearing our time resolution to 400 /is. No large 
samples were detected during the observation. We then sub- 
divided the pulse period into 256 phase bins for a time res- 
olution of 590 /is and examined the ^65000 single pulses. 
The effective rms was 70 mJy per phase bin, no single phase 
bin had a fiux density in excess of 5a. 



3.5 PSR J1617-5055 

The pulsar is very weak with a flux density of only 0.15 
mJy in our 3 hour observation. The pulsar is significantly 
scatter-broadened at this frequency and the peak flux is only 
1.5 mJy. We subdivided the pulse period into 128 phase bins 
and examined each of the ~ 150000 pulses. No single phase 
bin of any pulse was in excess of 400 mJy or 5x the noise 
level. This is consistent with Kaspi et al. (1998) who also 
reported a lack of giant pulses from this pulsar. 
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2.5 3 
LOG Flux Density (mJy) 



2.5 3 
LOG Flux Density (mJy) 



Figure 4. Cumulative probability distributions for two phase bins. The solid line is the expected distribution from Gaussian noise with 
an rms of 90 mJy. The right panel shows an extra solid line representing a power-law of slope -2.7. The dashed line is the reciprocal of 
the number of pulses. 



3.6 PSR B1706-44 

This pulsar was observed twice, each observation was 1.5 
hr in length. The integrated flux densities were 8.63 and 
7.70 mJy in the two observations. This pulsar has a pro- 
file which consists of a single component with a peak flux 
of ~100 mJy. Analysis of each 80^s sample revealed that 
many strong samples were detected in a small phase win- 
dow located far into the wings of the trailing edge of the 
profile. No such samples were detected on the leading edge 
and only a few at the peak of the pulse profile. It was clear 
that the large pulses were resolved on this timescale. We 
then divided the pulse period into 256 phase bins and the 
~46000 single pulses from each observation were examined. 
The rms per bin is 92 mJy. Figure 3 shows the pulse with 
the largest peak flux density. The flux density in this phase 
bin is 500 X the flux density in the integrated profile at this 
phase, but the integrated flux is only 4 times that of the in- 
tegrated profile (i.e. it would not be classified as a true giant 
pulse). The pulse is clearly resolved with a half width of ~1 
ms. This width must be intrinsic; the DM smearing is only 
120 ^,s and the scattering time at this freq uency is negligi- 
ble ( [Johnston, Nicastro fc Koribalski 1998| ). Of the ~93000 
pulses collected, 52 had a peak flux density in excess of 50 x 
the integrated flux density at a given phase. All these pulses 
are located between phases 0.06 and 0.095. 

Figure 4b shows the cumulative probability distribution 
of the intensities at phase 0.079. The mean flux density at 
this phase is 4 mJy, therefore only pulses with intensities 
greater than 80 x this can be detected at the 3-0" level. The 
distribution clearly deviates from that expected from gaus- 
sian noise and at high flux levels can best be described with 
a power law with index -2.7. This power law index is very 
similar to that seen in the Crab pulsar and PSR B1937-I-21. 
In contrast Figure 4a shows intensities for a bin on the rising 
edge of the proflle with the same integrated flux density as 
for phase 0.079. No large intensities are seen and the distri- 
bution is consistent with noise. 




Flux/Mean Flux 



Figure 5. Distribution of the integrated flux of all 96000 pulses 
of PSR B1706— 44 normalised by the mean flux (computed on a 
time-scale significantly shorter than the scintillation timescale) to 
remove the effects of interstellar scintillation. 



There are no true giant pulses seen in our data set 
for PSR B 1706-44. The maximum flux density of any sin- 
gle pulse was only ~4x the mean integrated flux density. 
Figure 5 shows the distribution of fluxes summing over all 
phase bins. The fluxes have been normalised to take into 
account the effects of interstellar scintillation. The width 
of the distribution is essentially determined by the receiver 
noise. However, there is an excess of counts at high fluxes 
- there are 782 pulses with flux densities greater than 2.4 x 
the mean flux density but only 490 pulses with flux densities 
less than -0.4 x the mean. This indicates the pulsar flux is 
likely log-normally distributed. 

It is our belief that the power-law distribution shown in 
Figure 3b is indicative of giant pulse behaviour. If this power 
law continues to larger fluxes, then 3.7 x 10^ rotations are 
needed before reaching a pulse with 20 times the mean flux 
density. This number of rotations is within a factor of 10 of 
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the giant pulse rate in both PSR B1937+21 and B1821-24 
but significantly higher than that of the Crab. 



4 DISCUSSION 

In the Vela pulsar, log-normal statistics are adequate to fit 
the distribution of flux across the bulk of the pulse profile 
(Kramer et al. 2001, Cairns et al. 2001) but the width of the 
distribution is significantly larger at the edge of the profile 
than in the middle. In addition to this, there are giant micro- 
pulses which occur well before the main pulse phase. These 
giant micro-pulses have a half-width of ~200 ^s, and are not 
at a fixed phase, but h ave an inherent 'jitter' of about 1 ms 
(Johnston et al. 2001). Their distribution is best described 
by a power- law (Kramer et al. 2001). We have found new 
examples of both these phenomena in our current survey. 

PSR B 1706-44 shows an additional example of giant 
micro-pulses. This time, however, the giant micro-pulses are 
located on the trailing edge of the pulse, and are somewhat 
wider then in Vela with a half-width close to 1 ms. Again 
there is some phase jitter as to the location of the pulse 
maximum. For this pulsar also, the distribution of fluxes is 
clearly power-law at high amplitudes. 

In PSR B1046-58 we clearly see large amplitude pulses 
on both the leading and trailing edges of the integrated pulse 
profile. The intensity distribution at these phases, however, 
is adequately described by a log-normal distribution with 
moderate width. The peak fluxes achieved are in excess of 
20 X the mean flux density but there is no evidence for a 
power-law tail to the distribution. In this regard they are 
similar to the fluctuations seen on the rising edge of the 
Vela pulsar. 

No giant micro- pulses are detectable in PSRs J 1105- 
6107, J1420-6048, B1509-58 or J1617-5055. However, the 
sensitivity to giant micro-pulses in these pulsars is not as 
good as for PSR B1706-44 or B1046-58. To demonstrate 
this, let's assume that any giant pulses would be similar 
to those seen in PSR B1706-44; i.e. they would have an 
intrinsic width of 0.5 ms, a power-law index of -3.0 and 
occur every 10^ rotations. The brightest giant pulse in the 
observing time for PSRs J1105-6107, J1420-6048, B1509-58 
and J1617-5055 would therefore be less than 500 mjy (in 
0.5 ms) and not detectable in the noise. Giant pulses could 
only be detected in these pulsars if they were intrinsically 
very narrow, and/or occured much more frequently than in 
either Vela or PSR B1706-44. 

We find, as in previous studies (e.g. Taylor et al. 1975), 
that the modulation of pulse intensity is larger in the wings 
of the profile of Vela-like pulsars. PSR B1046-58 is a good 
example of this behaviour: its phase resolved intensity dis- 
tributions are well described by a log-normal distribution 
whose width increases towards the proflle edge. The giant 
micro-pulses in Vela and now PSR B1706-44 appear to rep- 
resent a distinct pulse population, with emission phases well 
separated from the bulk of the integrated pulse profile. Both 
of these phenomena represent enhanced intensity fluctua- 
tions. 

Observations of the Crab pulsar at infra-red, optical 
and higher energies (e.g. Lundgren et al. 1995, Patt et al. 
1999, Romani et al. 2001) show that there are no detectable 
intensity fluctuations in the (incoherent) high energy emis- 



sion associated with the radio giant pulses. Thus we must 
conclude that giant pulses represent a variation in the co- 
herence of the particle distribution. An important question 
is why coherence fluctuations appear to be strongest at the 
pulse edges. 

One clue may be found in the connection of pulse width 
to altitude. Crudely, in a dipole polar cap of radius r, the 
width of the edge of the open zone maps as 5(j} oc r^^'^. If in- 
terpreted as the edge of a polar zone, the leading and trailing 
edges of the profiles where the modulation index increases 
represent locations 1.5-2x higher than that of the mean pul- 
sar emission. If we adopt the classical polar gap picture, with 
a pair formation front producing a dense plasma at a few 
stellar radii, then we can identify pulse intensity fluctuations 
with the growth of instabilities in the outflowing pair wind. 
Larger widths represent higher altitudes, and in this picture 
larger phase-space amplitudes for the instability-driven fluc- 
tuations, which have had more time to grow in the outflow. 
Alternative pictures exist, such as the slot-gap scenario of 
Arons (1983) for which the gap height is larger at the edge 
of the polar cap. Again there could be a plausible connection 
with increased instability. 

In both PSR B1706-44 and Vela the giant micro-pulses 
are located far from the peak of the pulse profile, repre- 
senting an emission altitude '-^4x higher. It is plausible that 
instability growth reaches a different regime at these heights. 
Clearly a true power-law intensity distribution suggests that 
the instabilities have left the (linear) stochastic growth 
regime (Cairns et al. 2001). Power-law behaviour suggests a 
causal connection between different intensity scales such as 
a cascade from one energy scale to another, or a self-similar 
coupling as in a saturated growth scheme. Such processes 
might develop late (higher) in the polar wind outflow. 

The giant micro-pulse emission is however well sepa- 
rated from the bulk of the radio emission in these objects, 
suggesting instead an independent origin. There is good ev- 
idence that for the Crab and possibly PSR B1821-24, the 
giant pulse radio emission is co-located with high energy 
emission in the outer magnetosphere, where strong power 
law X-ray components indicate sites of dense pair produc- 
tion. Vela and PSR B1706-44 are also 7-ray emitters; in 
the outer magnetosphere picture this high energy emission 
comes from the pole opposite to that viewed in the radio. 
They do not show strong X-ray pulses and so it is perhaps 
not surprising that we have not found giant pulse emission 
coincident with the 7-rays in these objects, as the Earth 
line-of-sight evidently does not sample regions of dense pair 
production. Nevertheless, there may be a high energy con- 
nection to the giant micro-pulse emission if the outer mag- 
netosphere above the radio pole is also actively producing 
pairs. The high energy emission from beyond the null charge 
surface would not be visible from this pole, but some pairs 
from the gap would be expected to flow inward past the null 
charge surface. We can speculate that this plasma mirrors 
in the converging field lines above the radio pole and that 
this mirrored, counter-streaming population would suffer in- 
stability growth and produce the giant micro-pulse compo- 
nents. 

A final question concerns the asymmetry of the giant 
micro-pulse component. Why does it lead in the case of 
Vela, but lag for PSR B 1706-44? For Vela, the radio pulse 
has a steep rise and slow fall suggesting the leading edge 
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of a cone; the giant micro-pulse component is at higher al- 
titudes on the same side. For PSR B1706— 44, one would 
then infer that the main radio pulse is the trailing edge of 
a cone, although this is morphologically less clear. If there 
is a high energy connection, then it is intruiging to note 
that the leading 7-ray pulse is stronger for Vela, while for 
PSR B1706— 44, the trailing component is stronger at most 
high energies. More examples are clearly needed to see if the 
dominance of leading versus trailing giant micro-pulses has 
a deterministic connection with other pulsar emission and 
if some global asymmetry in the magnetosphere geometry 
controls this choice. 



CONCLUSIONS 

We have found evidence for giant micro-pulses in PSR 

B1706-44 which arc very similar to those seen in the Vela 
pulsar. They are located in a small window of pulse phase, 
have a duration of ^1 ms and significant phase jitter. Their 
amplitude distribution is power-law and so may extend into 
true 'giant pulses' if observed for long enough. It is unclear 
whether the giant micro-pulses are just another manifesta- 
tion of pulsar 'weather' associated with the standard ra- 
dio emission from pulsars or whether they are more closely 
related to high-energy phenomena and the classical giant 
pulses. 

No giant pulses or any other giant micro-pulses were 
detected within our sensitivity limits in any of the other 
5 pulsars in our survey. Where it was possible to measure 
amplitude distributions, these were generally log-normal as 
surmised by Cairns et al. (2001). 
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